Maximal Gaussian entanglement achievable by feedback controlled dynamics 
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We determine a general upper bound for the steady-state entanglement achievable by continuous 
feedback for a system of any number of bosonic degrees of freedom. We apply such a bound to the 
specific case of parametric interactions - the most common practical way to generate entanglement in 
quantum optics - and point out optimal feedback strategies that achieve the maximal entanglement. 
We also consider the case of feedback schemes entirely restricted to local operations and compare their 
performance to the optimal, generally non-local, schemes. 

PACS numbers: 03.67.Bg, 02.30.Yy, 42.50.Dv 



The field of quantum control is central in the current 
rise of quantum technologies 0]. In particular, the 
control of the coherent resources of quantum states is 
an issue of major interest. Most valuable, and delicate, 
among such resources is certainly quantum entangle- 
ment, whose control is a primary requisite for quantum 
information and communication |3|-|6|]. This paper ad- 
dresses the question of how much entanglement can be 
generated by controlling the dynamics of a bosonic quan- 
tum system, and leads to the determination of optimal 
control schemes - achieving maximal entanglement - in 
relevant practical cases. In particular, we will consider 
systems subject to generic quadratic Hamiltonians and 
losses, and derive a bound on the maximal entanglement 
achievable, between specific bipartitions, by feedback 
schemes based on general continuous measurements 
and linear driving [7]. The class of dynamics and feed- 
back strategies covered in our study is very important in 
quantum optics, and is applicable to more general con- 
tinuous variable systems (ranging from atoms to nano- 
mechanical resonators). Being crucial for the implemen- 
tation of a number of quantum information protocols [8], 
the optimisation of the generation of continuous variable 
entanglement has been drawing considerable attention 
in recent years Since entanglement is not a linear 

figure of merit in the quantum state's parameters, one 
cannot tackle this optimisation with standard tools, like 
semi-definite programming j9fl, but rather requires the 
more detailed, specific analysis we shall present. 

Notation - We consider systems of N degrees of free- 
doms described by pairs of canonical operators: defin- 
ing a vector of operators x = {Q\,P\, —,^N/Pn) > ° ne nas 
= iCljt, where Q is the (2N)x(2N) symplectic form: 
= Vutl - (-iyi/2 - 5 jM1 [l + (-iy'1/2, in terms of 
Kronecker deltas 6^. Also, a - } = (<fy + ipj)/ sjl. 

For a system with such a phase-space structure we 
can define "Gaussian states" as the states with Gaussian 
Wigner functions. These states are completely deter- 
mined by the vector of means (x), and by the covariance 
matrix (CM) a, with entries cx* = ((AxjAxk) + (A£;AJjt)), 



where Ad = (d - (d)) for operator 6. The - always nec- 
essary - Robertson-Schrodinger uncertainty relation is 
also sufficient for Gaussian states to be physical [10]: 



a + z'Q > 



(1) 



We will consider Hamiltonians H that are at most of the 
second-order in x, so that their resulting free evolutions 
are affine in phase-space: H = (l/2)x T Hx - x T QBu(f), 
where the "Hamiltonian matrix" H is real and symmetric 
and B is real. The second term of H is a 'linear driving' 
proportional to a time-dependent input u(f): this term 
will describe the control exerted over the system. 

The system is considered to be open and such that 
each degree of freedom has its own channel to interact 
with the environment. Though thermal noise can also 
be treated along the lines we will present here, in this 
study we specialise for simplicity to pure losses, which 
are the main source of decoherence in quantum optical 
settings. We will thus assume a beam splitter-like ("ro- 
tating wave") interaction between each mode and the 
associated modes of the bath. Under the conditions set 
out above, the first moments of the canonical operators 
evolve according to d(x)/dt - A(x) + Bu(r), while the 
second moments obey 



da/df = Aa + oA T + 1. 



(2) 



Here, A = (QFJ - l)/2 is the "drift matrix", and 1 stands 
for the identity matrix with dimension clear from the 
context. We will only address stable systems, for which 
(A + A J ) < 0. Note that, for Gaussian states, these equa- 
tions describe the complete dynamics of the system. 

As customary in the context of feedback control, we 
will now assume that the degrees of freedom of the envi- 
ronment can be continuously monitored on time-scales 
which are short with respect to the system's response 
time fllll . The most general (efficient) measurement on 
the environment with outcomes continuous in time cor- 
responds to monitoring the operators (a T 1 + a f T), where 
the vector a = . . ,ajv) T contains all the annihila- 
tion operators of the system, and the complex matrix 
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T parametrises the measurement. These measurements 
(also known as 'general dyne' detections, see (ID) are 
very general, including heterodyne and homodyne de- 
tections as special cases, and define the broad setting of 
"continuous feedback" See for a description of 

the POVM giving rise to such measurements. In turn, Y 
defines the so called "unravelling matrix" U, given by 



l/l + Re[Y] Im[Y] 
2 I Im [Y] 1 - Re [Y] 



(3) 



The only conditions on Y are that U be symmetric and 
positive semi-definite. The outcome of the measure- 
ments on the environment is recorded as a "current" 
y = C<x) + where C = 2U 1/2 C and C jk = (6 2 /-u + 

5 2 (j-N),k)/ V2 for j, k e [1, . . . , 2N]. Finally, dw is a vector 
of real Wiener increments satisfying dwdw T = Idt (lj|. 
Clearly this treatment, like any feedback model, applies 
to systems where the output channels are open to exper- 
imental scrutiny like, e.g., light modes resonating in a 
cavity (where leaking light can be detected). The condi- 
tional evolution of the moments under such continuous 
measurements can be derived by standard techniques 
(Ito calculus). It amounts to a diffusive equation with a 
stochastic component for the first moments (x), and to 
a deterministic Riccati equation for the second moments 
fH. In our reasonings to follow, we will not make use of 
the details of such equations directly. We will be inter- 
ested in stable systems, and will determine the maximal 
entanglement achievable at steady state. Hence, all we 
need to remark is that a CM a is a stabilising solution Il2ll 
of the Riccati equation for the second moments if and 
only if d]: 



Aa + oA T + 1 > . 



(4) 



Together with Inequality Q), this relationship completely 
determines the set of stabilising solutions of our condi- 
tional dynamics. 

The final ingredient of the dynamics is the dependence 
of the linear driving u(f) on the history of the measure- 
ment record y(s) for s < t, which affects both first and sec- 
ond moments of the unconditional, 'average', evolution 
(whereas the second moments of the conditional states are 
unaffected by the linear driving), and closes the control 
loop. We will denote the unconditional state by q. Note 
that, for our class of dynamics, q is a statistical mixture of 
states with the same conditional CM a, obeying Inequal- 
ity (4j, and varying first moments. For Gaussian states, 
this implies that q can be obtained from a Gaussian state 
Qo with CM a and vanishing first moments by local op- 
erations and classical communication alone: q = L(qq), 
where L is some LOCC map. 

The typical aim of control over some time interval is 
to optimise the expected value of a cost function |TI. [T2|l . 
Our cost function will be the entanglement of Gaussian 
multi-mode steady states for bipartitions of 1 versus 



(N - 1) modes and 'bi-symmetric' bipartitions (i.e., in- 
variant under the permutation of local modes). Such an 
entanglement can be quantified by the logarithmic neg- 
ativity Ejv = — log 2 v_, where is the smallest eigen- 
value of (-aClaCi 1 ), Cl being the partial transposition of 
Q (TEH. Clearly, V- is not a quadratic cost function 
(i.e., it is not linear in a). This is why, albeit dealing with 
linear systems with Gaussian noise, we cannot resort 
to optimisation methods borrowed from classical LQG 
control problems OJ. 

General results - The main analytical result of this paper 
is presented here. Its proof may be found in appendix. 

Proposition 1 (Maximal entanglement) Let gbea steady 
state achievable by continuous Gaussian measurements and 
linear driving for a system of any number of bosonic modes 
subject to losses and to a Hamiltonian matrix H. The loga- 
rithmic negativity Em(q) of any 1 versus (N - 1) modes or 
bi-symmetric bipartition of q is bounded by: 



Em(q) < max 



1 



/-^ lo g 2 ( a H) 



(5) 



where {a -\ are the (strictly positive) eigenvalues of(-A - A T ) 
in increasing order, and A = j(QH - 1). 



Inequality 10 corresponds to 



v 2 _ > a.a 



(6) 



in terms of the smallest partially transposed symplectic 
eigenvalue of the Gaussian state Q. 

The bound above applies to both conditional and un- 
conditional states. In practice, only unconditional states 
are of interest since, although the first moments of the 
conditional states are in principle known, they fluctuate 
so fast (on the time-scale of the environment's dynamics) 
that the actual experimental state is the unconditional, 
average one. This is where the linear driving plays its 
crucial role in preserving the entanglement. Since the 
entanglement (for us, the logarithmic negativity) only 
depends on the second moments and decreases under 
LOCC, and since the second moments of the conditional 
states do not depend on the linear drive, the optimal 
choice for the linear driving is the one, always exist- 
ing, that keeps the first moments fixed (say, at zero). In 
this way, the linear drive's action guarantees that the 
unconditional state is at all times a conditional state - 
satisfying Inequality - with vanishing first moments. 
Hence, the optimal entangling strategy only depends on 
the optimal unravelling matrix U. 

Applications - Our theoretical result applies in general 
to all bosonic systems subject to losses and quadratic 
Hamiltonians. Here, we focus on optical modes os- 
cillating in a damped cavity and interacting through a 
parametric crystal or more general nonlinear media 
(a "non-degenerate, multi-frequency optical parametric 
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oscillator" [15]). Parametric interactions are the state of 
the art technology to generate continuous variable entan- 
glement. Also, optical bosonic systems can be interfaced 
with atomic systems Hill , so that the feedback scheme 
could be used to control atomic entanglement as well. 

The parametric interaction between modes j and k is 
described by the Hamiltonian xitfjpk + Pflk) We will 
assume equal interaction strengths % > between each 
pair of modes, consider a (n + ji)-mode bipartition, and 
describe analytically the scaling of the control of the en- 
tanglement with the number of modes n (we also define 
N = 2n). Our bound in this case is tight, and yields the 
actual optimal entanglement achievable by continuous 
filtering. Due to the symmetry of the system under the 
exchange of any two modes, the entanglement between 
the n-modes subsystems can be reduced to two-mode 
entanglement [18]: a local symplectic transformation ex- 
ists that turns the matrix A into an equivalent two-mode 
drift matrix A, plus a direct sum of irrelevant decoupled 
single-mode matrices The matrix A reads: 



( (« - ik 





nx 


> 







-(« - ik 





-nx 


1 


nx 





(n - lk 





~ 2 


I o 


-nx 





-(n - 1) X j 





For the system to be stable one must require: x < 2(N-i) 
(unstable systems, although in principle capable of gen- 
erating substantial entanglement, are in practice not con- 
trollable and certainly undesirable). As A is symmet- 
ric and invertible, the 'free' steady state CM oy can be 
promptly determined from Eq. Q: oy = -A~ 1 /2. Its log- 
arithmic negativity is given by i log 2 [(l + 2x)(l + 2(N - 
lk)]. Instead, the bound of Inequality 101 for any steady 
state CM a with continuous feedback control reads 



E N < -\ (log 2 (l - 2 X ) + log 2 [1 - 2(N - 1) X ]) 



(8) 



This upper bound is attained by the CM a ov t = 
_R T diag(a 2 / 1/«2/ l/«i/«i)K/ where R is the orthogonal 
transformation that diagonalises A and {a,} are the eigen- 
values of —2A in increasing order. This solution also satu- 
rates the Inequalities (2) and Q}. Both the free asymptotic 
entanglement and the optimal one under continuous fil- 
tering have thus been obtained analytically. Once the 
optimal achievable state is known as is the case here, 
the "optimal unravelling" U op t, and hence the optimal 
feedback scheme, can be straightforwardly derived since 
U op t = E(Aa opt + <JoptA J + 1)E T , where E = (2Ca opt - C) @]. 
For two modes, this rigorously proves that the schemes 
considered in Ref. [4] are indeed optimal. 
Local control - Such an optimal entanglement is in gen- 
eral achieved by filtering the system through global mea- 
surements on the environment, as no restrictions were 
assumed for the unravelling matrix U. This applies to sit- 
uations where the output channels of the two local sub- 
systems can be combined before being measured (like, 



e.g., for a parametric crystal in a cavity). We intend now 
to provide a lower bound on the entanglement achiev- 
able under local control, where the environmental de- 
grees of freedom pertaining to the separate subsystems 
cannot be combined, and compare it to the upper bound 
we obtained above. To this end, we will adopt direct 
(Markovian) feedback [7] and set u(t) - Fy(t). The un- 
conditional evolution of the system is then described by 

d<7/df = A' a + oA' J + D' , (9) 

with drift and diffusion matrices modified as A' = 
A + BFC and D' - 1 - C T F T B T - BFC + 2BFF J B J . We 
also choose a specific form of U and BF. Since in the 
free dynamics, governed by the drift matrix of Eq. (0, 
the quadratures p\ and pi are less noisy than c\\ and (fa, it 
is advantageous to monitor locally p\ and p2 and drive 
with the respective currents the quadratures ^2 and q\. 
However, due to the possible asymmetry of the two sub- 
systems for m + n, we have to consider different driving 
amplitudes \ii and \i% for their quadratures. All this 
corresponds to setting U33 = U44 = 1, V2(BF)24 = /J2, 
V2(BF)43 = and all other entries of LI and BF vanish- 
ing. We can then find the steady state solution of Eq. ^ 
as a function of the two feedback amplitudes and y.2, 
and evaluate its logarithmic negativity. It turns out that 
the maximum logarithmic negativity at steady state is 
attained for \i% — \i\njm. Hence, we are left with the en- 
tanglement depending on one parameter, over which we 
minimise numerically in the stable region, determined 
by (A' + A' J ) < 0. As a case of study, we have considered 
a system of 6 modes and summarised the results in Fig.Q] 
Because of the symmetry of the Hamiltonian, local con- 
trol is very close to optimal global control in the case of 
a balanced bipartition. However, the more unbalanced 
the bipartition, the more degraded the control, although 
numerics indicate that arbitrarily large entanglement can 
always be retrieved approaching the instability. 

Before concluding, let us further emphasise the useful- 
ness of feedback control by describing the practical case 
of two modes with interaction strength to loss factor 
ratio x — 0.45. Without control, this system would gen- 
erate 0.93 ebits of logarithmic negativity at steady state. 
The optimal feedback control would rise this value to 
3.32 ebits. The Markovian local control discussed here, 
instead, allows one to reach 2.12 ebits: a remarkable im- 
provement over the case with no control showing that, 
in this instance, about half of the entanglement retriev- 
able by global measurements can be recovered from the 
environment by local measurements. 
Conclusion - We derived a bound on the entanglement 
achievable, at steady state and for various bipartitions, 
in multimode linear bosonic systems under continuous 
feedback control. When applied to the practical case 
of symmetric parametric interactions, our bound also 
allows one to determine the measurement strategy max- 
imising the steady-state entanglement, which is relevant 
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FIG. 1: Squared symplectic eigenvalue vz at steady state for 
a system of 6 modes (P_ — » implies infinite entanglement). 
Green (lighter) curves depict nvr_ in the absence of control (from 
top to bottom: 1:5, 2:4, and 3:3 modes bipartition); bleu (darker) 
curves refer to numerically optimised local feedback (from top 
to bottom: 1:5, 2:4, and 3:3 modes bipartition); the red curve is 
the analytical lower bound (|6). 

to optimise the experimental generation of continuous 
variable entanglement, and hence useful for countless 
quantum information protocols [8]. More generally, our 
investigation yields a technique for the optimisation of 
nonlinear figure of merits in bosonic quantum systems, 
with a broad range of applications in quantum informa- 
tion processing and state engineering. 
Appendix - Proof of Proposition 1. Henceforth, \v) will 
stand for a unit vector in the phase space T and { A^} 

({A 1 }) will be the 2N increasingly-ordered (decreasingly- 
ordered) eigenvalues of an N-mode CM a. For each \v), 
one can define the unit vector \w) = Cla l ^ 2 \v)l ^l(v\a\v), 
such that (v\o 1/2 \w) = ( since Q = -Q T ) and 

v 2 _ > min {v\a\v){w\a\w) = A\a\ . (10) 

with the min taken over \v), \w) satisfying {v\a ll2 \w) = 0. 

We will further denote by \vi) the eigenvectors corre- 
sponding to the increasingly ordered eigenvalues of a: 
a\vj) = Aj\vj). Then, by using the Robertson Schrodinger 
Inequality and the 'Poincare Inequality' fl9ll , one can 
show that a vector \w) must exist in Or; c (defined as the 
subspace spanned by the k orthogonal vectors Q|c/ C )) for 
which (w\o\w) < Aj[, and such that 

A\a[ > 1. (11) 

Now, let Ock, be a conditional CM at steady state ob- 
tained under continuous measurements, pure losses and 
a Hamiltonian matrix H. Applying Inequality (UJl to the 
eigenvectors corresponding to A] and A^, one has for the 
two largest eigenvalues Aj and A^ of a: 

A\A l 2 < -L , (12) 
a[a 2 



where {a 1 } are the (strictly positive) eigenvalues of (-A - 

A T ) in increasing order. The chain of Inequalities (10) . 
(TTjl and {12)1 leads to I0 for the partially transposed 
symplectic eigenvalue of the conditional state. 

Finally, as we have seen previously, g = L(qq), where 
L is a LOCC operation and £0 a Gaussian state with a 
CM which is a stabilising solution of 10. Hence Eyv(f) = 

Ejv(L(p )) < Ejv(<?o) ^ max [0, - log 2 (flJ«2)/2]/ where I©, 
the formula Ey = - log 2 (v_), and the monotonicity of Ejv 
under LOCC |2(J have been invoked. □ 
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